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We study superconductivity in the two-dimensional 16-band d-p model extracted from a tight- 
binding fit to the band structure of LaFeAsO, using the random phase approximation. When the 
intraorbital repulsion U is larger than the interorbital one U' , an extended s-wave (s±-wave) pairing 
with sign reversal of order parameter is mediated by antiferromagnetic spin fluctuations, while when 
U < U' another kind of s-wave (s++-wave) pairing without sign reversal is mediated by ferro-orbital 
fluctuations. The s++-wave pairing is enhanced due to the electron-phonon coupling and then can 
be expanded over the realistic parameter region with U > U' . 



I. INTRODUCTION 

The recently discovered iron-based superconductors 1 ^ 
RFePnOi.^Fz (R=Rare Earth, Pn— As, P) with a 
transition temperature T c exceeding 50K— ~— have at- 
tracted much attention. The F nondoped compound 
LaFeAsO exhibits the structural transition from tetrago- 
nal (P4/nmm) to orthorhombic (Cmma) phase at a tran- 
sition temperature T =155K and stripe-type antiferro- 
magnetic order at T — 134K with a magnetic moment 
~ 0.36^5 8 at low temperature. With increasing F dop- 
ing, the system becomes metallic and the antiferromag- 
netic order disappears 2 , and then, the superconductivity 
emerges for x ~ 0.1 with T c ~ 26K. Rare-earth substitu- 
tion compounds exhibit superconducting transition with 
higher T^-~—. Specific features of the systems are two- 
dimensionality of the conducting Fe2As2 plane and the 
orbital degrees of freedom in Fe 2+ (3d 6 )^. The pairing 
symmetry together with the mechanism of the supercon- 
ductivity is one of the most significant issues. 

The NMR Knight shift measurements revealed that 
the superconductivity of the systems is the spin-singlet 
pairing 9 ' 10 . Fully gapped superconducting states have 
been predicted by various experiments such as the pen- 
etration depthii, the specific heati£, the angle resolved 
photoemission spectroscopy (ARPES)i2r— and the im- 
purity effect on T ^ 10 ^ 6 . In contrast to the above men- 
tioned experiments, the NMR relaxation rate shows the 
power low behavior \jT\ oc T 3 below T c 17 , suggesting 
the nodal or highly anisotropic gap structure. The other 
NMR measurements— however, revealed \jT\ cx T 6 be- 
low T c and there is still controversy. 

Theoretically, the first principle calculations have pre- 
dicted that the nondoped system is metallic with two or 
three concentric hole Fermi surfaces around the T point 
(k = (0,0)) and two elliptical electron Fermi surfaces 
around the M point (k = (w,tt))^~— . Mazin et al. 
suggested that the spin-singlet extended s-wave pairing 
whose order parameter changes its sign between the hole 
pockets and the electron pockets (s±-wave) is favored due 
to the antiferromagnetic spin fluctuations 2 ^. According 
to the weak coupling approaches based on multi-orbital 



Hubbard models 2 ^—, the s±-wave pairing or the d xy - 
wave pairing is expected to emerge. It is shown that 
the s±-wave pairing is realized also in the strong cou- 
pling region by the mean field study based on the t-Ji-J% 
model— and the exact diagonalization study based on 
the one-dimensional two-band Hubbard model 3 ^. 

Generally speaking, the details of the band structure 
and the Fermi surface are crucial for determining the 
pairing symmetry. In the 5-band Hubbard model origi- 
nally introduced by Kuroki et al., the energy bands ob- 
tained by reproduce those obtained by the density func- 
tional calculation very well 2 ^. In this model, however, the 
spatial extensions of the Fe 3d like Wannier orbitals are 
different from each other— and the resulting intra-orbital 
terms of the on-site Coulomb interaction are strongly or- 
bital dependent 38 . In addition, since the model explic- 
itly includes the transfer integrals up to the fifth nearest 
neighbor sites 2 ^, one should take the off-site Coulomb 
interaction, which is considered to be about 0.5eV be- 
tween the nearest neighbor sites, into account to ensure 
the consistency of the model 3 ^.. On the other hand, in 
the effective model which includes both the Fe 3d orbitals 
and the As Ap orbitals, so called d-p model, the spatial ex- 
tensions and the differences of those between the orbitals 
will be considerably reduced 3 - 7 -. Due to these facts, in 
the d-p model, it is expected that the intra-orbital terms 
of the on-site Coulomb interaction for each orbitals have 
almost the same values and the off-site Coulomb interac- 
tion are negligible. Therefore, theoretical studies based 
on the d-p model, are highly desired. 

In the previous papers 3 ^—, we have investigated the 
electronic states of the Fe2As2 plane in iron-based super- 
conductors on the basis of the two-dimensional 16-band 
d-p model which includes the Coulomb interaction on a 
Fe site: the intra- and inter-orbital direct terms U and 
U', the Hund's coupling J and the pair-transfer J'. Using 
the random phase approximation (RPA), we have found 
that, for a larger value of J, the most favorable pairing 
symmetry is s±-wave, while, for a smaller value of J, it 
is d X y-w&ve. 

The present paper is a full paper to our previous 
papers 3 ^— with some numerical improvements 41 . In the 
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present paper, we investigate the superconductivity in 
the wider parameter space by treating U,U', J and J' as 
independent parameters in contrast to the previous study 
under the condition that U = U' + 2 J and J = J' based 
on the two-dimensional 16-band d-p model. Solving the 
superconducting gap equation with the pairing interac- 
tion obtained by using the RPA, we find that two kinds 
of the s-wave superconducting states appear. As above 
mentioned, the s±-wave superconducting state emerges 
near the incommensurate spin density wave (ISDW) with 
q ~ (n, tt) phase. In addition, for U < U', the s-wave su- 
perconducting state appears around the ferro-orbital or- 
dered phase. The order parameter for this s-wave state 
dose not change its sign in k space. We refer to this 
s-wave state as the s++-wave state, hereafter. 



II. MODEL AND FORMULATION 

First of all, we perform the density functional calcu- 
lation for LaFeAsO with the generalized gradient ap- 
proximation of Perdew, Burke and Ernzerho^ by us- 
ing the WIEN2k package 43 , where the lattice parameters 
(a = 4.03268A, c = 8.7411lA) and the internal coordi- 
nates (zLa = 0.14134, zas — 0.65166) are experimentally 
determined 4 ^. The crystal structure of Fe2As2 layer is 
shown in Fig. [T](a). Since As atoms are tetrahedrally 
arranged around a Fe atom, there are two distinct Fe and 
As sites in the crystallographic unit cell (see Figs. [U(a), 
(b)). Considering these facts, we then derive the two- 
dimensional 16-band d-p mode l 39 ' 40 , where 3d orbitals 
(e? 3z 2_ r 2, d x 2_ y 2, d xy , d yz , d zx ) of two Fe atoms (Fe 1 =A, 
Fc 2 =B) and Ap orbitals (p x , p y , p z ) of two As atoms are 
explicitly included. We note that x, y axes are directed 
along second nearest Fe-Fe bonds (see Fig. Q](b)). 

The total Hamiltonian of the d-p model is given by 
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H = H + H h 



(1) 



where Hq and H m t are the nonintcracting and interact- 
ing parts of the Hamiltonian, respectively. The nonin- 
teracting part of the d-p model is given by the following 
tight-binding Hamiltonian, 
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where dn a is the annihilation operator for Fe-3d elec- 
trons with spin a in the orbital £ at the site i and Pi m a 
is the annihilation operator for As-4p electrons with spin 
a in the orbital m at the site i. In eq. @, the transfer 
integrals t dd lt , t pp jmm ,, t dp j/ m and the atomic ener- 
gies e d , e v m are determined so as to fit both the energy 
and the weights of orbitals for each band obtained from 
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FIG. 1. (Color online) Crystal structure of Fe2As2 layer . 
Small and large balls represent Fe and As atoms, respectively. 
The solid line represents the unit cell. It is noted that As 
and As 2 denote the As atoms on the upper side and on the 
lower side of the Fe2As2 layer, respectively. 



the tight-binding approximation to those from the den- 
sity functional calculation 4 ^. Similar models have been 
used by the other author a 34 ' 46 ' 47 but the model param- 
eters are different from ours. The doping concentration 
x corresponds to the number of electrons per unit cell 
n = 24 + 2x in the present model. 

Now we consider the effect of the Coulomb interaction 
on Fe site. The interacting part of the Hamiltonian is 
given as follows, 

ff i»t = ^EEE 4 ery di es d U9 d U a 



f d f d- 
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where U and U' are the intra- and inter-orbital direct 
terms, respectively, and J and J' are the Hund's cou- 
pling and the pair-transfer, respectively. For the isolated 
atoms, the relations between Coulomb matrix elements 
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U = U' + 2 J and J — J' are derived due to the rotational 
invariance of the Coulomb interaction and the reality of 
the wave functions, respectively 4 ^. For the atoms in the 
crystal, however, the relation is not satisfied generally 
due to the crystallographic effects and the many body 
effects due to the Coulomb interaction and the electron- 
phonon coupling which will be discussed later. Therefore, 
we treat U, U' , J and J' as independent parameters in 
the present paper. 

Within the RPA 4 ^— , the spin susceptibility x s (q) an d 
the charge-orbital susceptibility x c (q) are given in the 
50 x 50 matrix representation as follow s 39 ' 40 , 



on-site energy 



X s (q) = (l-X W (q)5)- 1 x W (q), 
X c (q) = (i + X ( ° ) (q)C)- 1 x (0) (q) 
with the noninteracting susceptibility 
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where /x, ^ (=1-16) are band indexes, a, /3 (=A, B) rep- 
resent two Fe sites, I represents Fe 3d orbitals, (k) is 
the eigenvector which diagonalizes H eq. ([5} , „ is the 
corresponding eigenenergy of band \x with wave vector k 
and /(e) is the Fermi distribution function. In eqs. (4) 
and (5), the interaction matrix 5* (C) is given by 
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In the weak coupling regime, the superconducting gap 
equation is given b y 39 i 40 
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where A^, (k) is the gap function and V^f 2 1 ^ (q) is 
the effective pairing interaction 52 . Within the RPA4 9 - 
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TABLE I. Tight-binding parameters (in units of eV) for the 
d-p Hamiltonian eq. ©. It is noted that we define the d-p 
hopping and the in-plane p-p hopping parameters along x- 
axis. 



compounds^. For simplicity, we set x = 0.1 and T = 
0.02eV in the present study. We use 32 x 32 k points in 
the numerical calculations for eqs. (H])-©, and also use 
the fast Fourier transformation (FFT) to solve the gap 
equation eq. ([5]). Here and hereafter, we measure the 
energy in units of eV. 
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where r) = § for the spin-singlet state and T] = — \ for 
the spin-triplet state. The gap equation flS) is solved 
to obtain the gap function A"^(k) with the eigenvalue 
A. At T = T c , the largest eigenvalue A becomes unity. 
In the present paper, we only focus on the case with 
x = 0.1, where the superconductivity is observed in the 



III. CALCULATED RESULTS 

A. Band Structure 

We show the band structure obtained from the d- 
p tight-binding Hamiltonian eq. ©, where the tight- 
binding parameters are listed in Table IT] 15 , together with 
that obtained from the density functional calculation in 
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FIG. 2. (Color online) (a) The band structure obtained from 
the d-p model eq. ([2} (solid line) and that obtained from the 
density functional calculation (cross) for x = 0.1. (b) Fermi 
surface obtained from the d-p model for x — 0.1. The solid 
and dashed lines show the Fermi surfaces which have mainly 
dy Z ,d zx and d x 2_ y 2 orbital character, respectively, (c) The 
DOS obtained from the d-p model for x = 0.1. Upper panel: 
total DOS, Middle panel: partial DOS of the Fe 3d orbitals, 
Lower panel: partial DOS of the As 4p orbitals. The inset 
of the middle panel shows the DOS near the Fermi level. We 
note that the Fermi level set to on the energy axis. 



the Fig. [5](a). The result of our density functional calcu- 
lation is similar to that previously reported by the other 
authors ^ 23 ' 25 . It is found that the former reproduces 
the latter very well. We note that the weights of orbitals 
also agree very well with each other (not shown). Due 
to the weak crystalline electric field from the As 3 ~ ions 
tetrahedrally arranged around a Fe atom and the strong 
hybridization between the Fe 3d orbitals, the resulting 



energy bands have very complicated structure. 

The Fermi surface for the d-p tight-binding Hamilto- 
nian is shown in Fig. [H(b), where we can see nearly circu- 
lar hole pockets around the T point and elliptical electron 
pockets around the M point. These results are consistent 
with the previous first principle calculations^^— . 

The density of states (DOS) obtained by the d-p tight- 
binding Hamiltonian eq. © is shown in Fig. [5](c). It 
is found that the dominant contribution near the Fermi 
level comes from Fe 3d orbitals and the contribution of 
As 4p orbitals is small but is not negligible. We show the 
partial DOS of Fe 3d orbitals and that of As Ap orbitals 
in the middle panel and the lower panel of Fig. [2](c), re- 
spectively. The d yz , d zx and d x 2_ y 2 states comprise the 
large part of the DOS near the Fermi level, while, the 
c?3 Z 2_ r 2, d xy states occupy the small one and are com- 
parable with the p x , p y and p z states. The d yz and d zx 
states at the Fermi level are larger than the d x 2_ y 2 ones 
and this corresponds to the fact that the electron pock- 
ets have d yz , d zx and d x 2_ y 2 orbital characters, while, 
the hole pockets have only d yz and d zx orbital charac- 
ters. However, the d x 2_ y 2 states have large values just 
below the Fermi level as shown in the inset of the mid- 
dle panel of Fig. [2] (c). This is due to the hole band 
near the T-point just below the Fermi level. Therefore, it 
is anticipated that the d yz , d zx and d x 2_ y 2 orbitals play 
significant roles to determine the magnetic, orbital and 
superconducting properties. 



B. RPA Results for U > U' 

In this subsection, we concentrate our attention on the 
case with U > U'. We set the typical parameters as U = 
1.71, U' = 1.4 and J = J' = 0.1, where the condition 
for the superconducting transition A = 1 is satisfied as 
mentioned below (see Fig. [3] (d)). 

The several components of the spin susceptibility x s (q) 
given in eq. (HJ are plotted in Fig. [3] (a) . The spin sus- 
ceptibility is enhanced due to the effect of the Coulomb 
interaction. It is found that the most dominant compo- 
nent is the d x 2_ y 2 diagonal component and the incom- 
mensurate peaks around the M point are observed as 
reflecting the nesting between the hole pockets and the 
electron pockets. As mentioned before, the hole band 
which has mainly d x 2_ y 2 orbital character exists just be- 
low (~ O.OleV) the Fermi level and contributes to the 
large value of the DOS (see Fig. [2] (c). Therefore, the 
d x 2_ y 2 diagonal component of X s (q) becomes most dom- 
inant at finite temperature T = 0.02eV (> O.OleV). The 
result is consistent with the RPA results based on the 
5-band Hubbard model 2 ^. 

The several components of the charge-orbital suscepti- 
bility x c (q) given in eq. ([5]) are plotted in Fig. [3](b). In 
contrast to the case with the spin susceptibility, the off- 
diagonal component oid x 2_ y 2 —d yz which corresponds to 
the transverse orbital susceptibility becomes most domi- 
nant and shows peaks around the M point together with 
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FIG. 3. (Color online) (a) Several components of the spin susceptibility x s (q), (b) the charge-orbital susceptibility X c (ci) an d 
(c) the effective pairing interaction V'(q) in q space for U = 1.71, U' = 1.4 and J — J' = O.leV. It is noted that we number 
the orbitals as follows: d 3z 2_,,2(l), d x 2_ y 2(2), d xy (3), d yz (4), d zx (5). (d) (7 dependence of the eigenvalues of the gap equation 
for several symmetries. 
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FIG. 4. (Color online) (a) Several components of the spin susceptibility x" (q), (b) the charge-orbital susceptibility x c (c[) an d 
(c) the effective pairing interaction t^(q) in q space for U = 0.4, U' = 1.15 and J — J' — O.leV. It is noted that (d) U' 
dependence of the eigenvalues of the gap equation for several symmetries. It is noted that the legends are the same as Fig. [3] 



those at the T point. It is noted that for U > U' the 
spin fluctuations dominate over the charge-orbital fluc- 
tuations as shown in Figs. [3] (a) and (b). 

The several components of the effective pairing inter- 
action V(q) for the spin-singlet state given in eq. © are 
plotted in Fig. [3] (c). Since the largest eigenvalue A is 
always spin-singlet state in the present study, we show 
the effective pairing interaction only for the spin-singlet 
state. Since in the case for U = 1.71, U' = 1.4 and 
J = J' = 0.1, the spin fluctuations dominate over the 
orbital fluctuations as mentioned above, the structures 
of V(q) are similar to those of the spin susceptibility. 

Substituting V(q) into the gap equation eq. ©, we 
obtain the gap function A(k) with the eigenvalue A. In 
Fig. El (d), the eigenvalues A for various pairing symme- 
tries are plotted as functions of U for fixed values of U', 
J, J' . With increasing U, A monotonically increases and 
finally becomes unity at a critical value U c above which 
the superconducting state is realized. For U' — 1.4 and 
J = J' = 0.1 the largest eigenvalue A is for the s-wave 
symmetry and U c = 1.71. The second largest eigenvalue 
is for d xy -wave symmetry and the eigenvalue for the d xy - 



wave symmetry increases as J increases for U > U'. 



C. RPA Results for U < U' 

In this subsection, we concentrate our attention on the 
case with U < U'. We set the typical parameters as U = 
0.4, U' = 1.15 and J = J' = 0.1, where the condition 
for the superconducting transition A = 1 is satisfied as 
mentioned below (see Fig. |U(d)). 

The several components of the spin susceptibility x s (q) 
given in eq. Q are plotted in Fig. [4] (a). In contrast to 
the case with U > U' (see Fig. [3J (a)), the off-diagonal 
element d x 2_ y 2 —d yz is most dominant owing to the inter- 
orbital direct term U' > U . 

The several components of the charge-orbital suscep- 
tibility x c (q) given in eq. §5§ are plotted in Fig. 2] (b). 
In contrast to the case with the spin susceptibility, the 
diagonal component of d x 2_ y 2 becomes most dominant 
and shows peaks around the T point. It is noted that 
for U < U' the charge-orbital fluctuations, which corre- 
sponds to the fluctuations near the ferro-orbital ordered 
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state realized in the large U' regime as mentioned later 
(see Fig. [5]), dominate over the spin fluctuations as shown 
in Figs. H (a) and (b). 

The several components of the effective pairing inter- 
action V"(q) for the spin-singlet state given in eq. ([9]) are 
plotted in Fig. g] (c). Since for U = 0.4, U' = 1.15 and 
J = J' = 0.1, the charge-orbital fluctuations are larger 
than the spin fluctuations, the diagonal components of 
V(q) are always negative in q space. 

Substituting V(q) into the gap equation eq. <JSj> , we 
obtain the gap function A(k) with the eigenvalue A. In 
Fig. S](d), the eigenvalues A for various pairing symme- 
tries are plotted as functions of U' for fixed values of U, 
J and J'. With increasing U', A monotonically increases 
and finally becomes unity at a critical value U' c = 1.15 
above which the superconducting state is realized. Sim- 
ilar to the case of U > U', the largest eigenvalue A is 
for the s-wave symmetry but the superconducting gap 
structure is significantly different from that for U > U' 
as shown below. 



D. Gap Functions 

First, we discuss the gap functions in the case with 
U > U'. Fig. [5] shows the diagonal components of the 
gap function A(k) for U = 1.71, U' = 1.4 J = J 1 = 0.1. 
Figs. [5](a)-(d) show the gap functions in the orbital rep- 
resentation and Figs. [5] (e)-(h) show those in the band 
representation. We note that the energy bands are num- 
bered as descending energy. It is found that the gap 
function has the s-wave symmetry and the most domi- 
nant component is the d x 2_ y 2 diagonal component. We 
find that the gap functions in the band representation 
have different signs between the electron pockets and the 
hole pockets without any nodes on the Fermi surfaces 
(s ± -wave symmetry It is noted that the diagonal 
components of the gap function in the orbital represen- 
tation, except for the d xy component, also change those 
signs in k space. The absolute values of the gap functions 
on the Fermi surfaces are almost isotropic but largely de- 
pend on the energy bands; those on the electron pockets 
of the 13th and 14th bands are twice or more larger than 
those on the hole pockets of the 11th and 12th bands. 
This is because the d x 2_ y 2 component, which has dom- 
inant contribution in % s (q) as shown in Fig. [3] (a) , for 
the 13th and 14th bands is larger than that for the 11th 
and 12th bands. We note that the 10th band (hole band) 
with the largest d x 2_ y 2 component has the largest abso- 
lute value of the gap function, although the Fermi level 
is just above the 10th band and does not cross it for 
x = 0.1. 

Next, we discuss the gap functions in the case with 
U < U'. Fig. |5]shows the diagonal components of the gap 
function A(k) for U = 0.4, U' = 1.15 J = J' — 0.1. Figs. 
|6](a)-(d) shows the gap functions in the orbital represen- 
tation and Figs. IH(e)-(h) show those in the band repre- 
sentation. The diagonal components of the gap function 
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FIG. 5. (Color online) (a), (b), (c), (d) The diagonal compo- 
nents of the gap function A(k) in the orbital representation 
and (e), (f), (g), (h) those in the band representation for 
U = 1.71, U' = 1.4 and J = J' = O.leV. The solid and 
dashed lines represent the Fermi surfaces and the nodes of 
the gap function, respectively. 



in the orbital representation have no sign change in the 
k space due to the diagonal components of V'(q) < as 
shown in Fig. [4] (c) . We call this s-wave state as the s++- 
wave state. The gap function in the band representation, 
however, has sign change between the Fermi surfaces and 
line nodes on the 14th band Fermi surface. These facts 
reflect that the sign change of the gap function in the 
orbital representation between the d x 2_ y 2 diagonal com- 
ponent and the other orbital diagonal components. The 
11th band and 12th band Fermi surface has mainly d yz 
and d zx orbital character, while, the 13th band Fermi 
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FIG. 6. (Color online) (a), (b), (c), (d) The diagonal compo- 
nents of the gap function A(k) in the orbital representation 
and (e), (f), (g), (h) those in the band representation for 
U = 0.4, U' = 1.15 and J = J' = O.leV. The solid and 
dashed lines represent the Fermi surfaces and the nodes of 
the gap function, respectively. 
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FIG. 7. (Color online) (a), (b), (c), (d) The diagonal compo- 
nents of the gap function A(k) in the orbital representation 
and (e), (f), (g), (h) those in the band representation for 
U = 0.4, U' = 1.18 and J = 0.1, J' = OeV. The solid and 
dashed lines represent the Fermi surfaces and the nodes of the 
gap function, respectively. 



surface has mainly d x 2_ y 2 orbital character. Therefore, 
the gap function has different sign between the hole pock- 
ets and the 13th band electron pocket. The 14th band 
electron pocket has mainly d yz and d zx orbital character 
away from the Brillouin zone boundary, while d x i^ y i or- 
bital character on the 14th band electron pocket is com- 
parable with d yz and d zx one near the Brillouin zone 
boundary. Thus, the gap function on the 14th band elec- 
tron pockets has plus sign near the zone boundary and 
minus sign away from the zone boundary. By the sim- 
ple mean field analysis of the pair transfer term of the 



interacting part of the Hamiltonian eq ([3]), 

It is shown that the pair transfer J' > favors the sign 
change between the diagonal components of the gap func- 
tion in the orbital representation. 

In fact, we also examine the case with J' = and we 
find that the s+_|_-wave state without sing change between 
the orbitals is realized for U < U'. We show the gap 
function for U = 0.4, U' = 1.18, J = 0.1, J' = in Fig. 
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[7J It is found that the all diagonal components of the gap 
function in the orbital representation have the same sign 
and those in the band representation have no sign change 
between all Fermi surfaces. Therefore, it is considered 
that the sign change of the gap function between the the 
d x 2_ y 2 diagonal component and the others is due to the 
pair transfer J /53 ' 54 . 

It is helpful for understanding the difference between 
the s±-wave state and the s++-wave state in more detail 
to consider the gap function in the real space. For s±- 
wave state, the on-site pairing is comparable with the 
nearest neighbor and/or the next nearest neighbor one. 
On the other hand, for the s++-wave state, the on-site 
pairing is dominant and the off-site pairings are negligibly 
small as compared to the on-site pairing. 

Here we discuss the reason why the on-site part of the 
gap function for s±-wave state is large (especially in the 
d x 2_ y 2 diagonal component) even though the most dom- 
inant component of the effective interaction is always re- 
pulsive in q space (see Fig. (c)), and then the on-site 
effective interaction is repulsive. When we perform the 
Fourier transformation of the gap equation eq. (|8j). the 
on-site part of the left hand side is proportional to the 
on-site gap function, while that of the right hand side 
is given by the product of the on-site effective interac- 
tion (x (— 1)) and the on-site anomalous Green's function 
which is proportional to the q summation of the gap- 
function times the single-particle spectral weight times 
the thermal factor. In the case with the d x 2_ y 2 diagonal 
component, the on-site gap function is negative as the 
negative contribution of the gap function in q space is 
much larger than the positive one as shown in Fig. [5] 
(a). On the other hand, the on-site anomalous Green's 
function becomes positive as the single-particle spectral 
weight of the d x 2 _ y 2 hole band is very large around the V 
point where the gap function is positive as compared to 
that of the electron band around the M-point where the 
gap function is negative. Then, the gap equation can be 
satisfied with the large value of the on-site gap function 
against the repulsive on-site effective interaction. 

When the doping x increases, the Fermi level rises 
apart from the d x 2_ y 2 hole band, and then the effect of 
the hole band decreases resulting in the decrease in the 
on-site gap function as well as the decrease in the super- 
conducting transition temperature (not shown). Such 
doping dependence of the on-site gap function has re- 
cently been observed in the 5-band Hubbard model 56 . 
On the contrary, in the s++-wave state, the on-site pair- 
ing is always dominant almost independent of the doping 



E. Phase Diagram 

The phase diagram on U'-U plane for J — J' = 0.1 
is shown in Fig. El where the magnetic and charge- 
orbital instability is determined by det(l — x'°^(q)5) = 
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FIG. 8. (Color online) The phase diagram on U'-U plane for 
J = J' = O.leV at x = 0.1, T = 0.02eV. The solid and dotted 
lines show the s++-wave and the s±-wave superconducting 
instabilities, respectively. The dashed and dot-dashed lines 
show instabilities towards the incommensurate spin density 
wave and the ferro-orbital order, respectively. 



conducting instability is determined by A = 1 as men- 
tioned before. The ISDW with q ~ (tt, it) appears in 
the large U region, while, the ferro-orbital order ap- 
pears for U < U— [see also Figs. [3] (a), (b) and Figs 
HI (a), (b)]. It is noted that on the phase boundary where 
the charge-orbital instability takes place, the longitudi- 
nal orbital susceptibility (^(l^eeew diverges, while, the 

charge susceptibility Y,e,e l ,a,p(x c ( c L))ili>e' dose not . The 
s±-wave pairing is realized near the ISDW due to the spin 
fluctuations, while, the s++-wave pairing is realized near 
the ferro-orbital ordered phase due to the charge-orbital 
fluctuations, where we regard the superconducting states 
as the s ++ -wave states if d yz , d zx and d x 2_ y 2 diagonal 
components of the gap function have no sign change in 
k space and as the s±-wave states if not. The way to 
determine whether the superconducting state is the s±- 
wave state or the s+_|_-wave state is not unique. This is 
because the s±-wave and the s++-wavc state arc same 
symmetry (Ai s ) and the change between s ++ -wave and 
the s±-wave state is crossover. In fact, as U increases, 
the on-site paring decreases, while, the off-site pairing 
increases continuously. At U ~ 1.25, the nodes appear 
around the M-point for the d x 2_ y 2 diagonal component 
and those approaches the T-point as U increases. It is 
noted that we also obtain the phase diagram on U'-U 
plane for J = J' = 0.25 and find that the phase diagram 
is essentially the same as that for J = J' = 0.1 except 
that the magnetic and the s±-wave superconducting in- 
stabilities are slightly enhanced by the larger value of the 
Hund's coupling J. 
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F. Effects of Electron-Phonon Coupling 

In this subsection, we discuss the effects of the electron- 
phonon coupling. By performing the group theoreti- 
cal analysis for LaFeAsO, it is found that there are 
14 kinds of the optical phonon modes at the T point: 
2Ai g +2Big+4E g +3A2 U +3E u . Here, we concentrate on 
the A\ g mode in which La and As ions oscillate along 
the c-axis. The A\ g phonon dose not break the sym- 
metry of the orbital and the resulting electron-phonon 
coupling matrix g is diagonal in the orbital representa- 
tion. Within the RPA, the charge-orbital susceptibility 
X c (q) including the effects of both the electron-electron 
and the electron-phonon coupling is obtained by replac- 
ing U with U - 2U ph and 2U' - J with 2(U' - U ph ) - J 
in eqs. §5§ and (J7]), where U ph — 2g 2 /uj Alg , uj Alg is the 
frequency of the A\ g phonon and we neglect the orbital- 
and q-dependence of the electron phonon interaction. It 
is found that the inter-orbital direct term U' which en- 
hances the orbital fluctuations is harder to be reduced 
by the electron-phonon coupling than the intra-orbital 
direct term U . As a result, the orbital fluctuations are 
relatively enhanced by the electron-phonon coupling as 
compared to the spin fluctuations. 

IV. SUMMARY AND DISCUSSION 

In summary, we have investigated the pairing symme- 
try of the two-dimensional 16-band d-p model by using 
the RPA and have obtained the phase diagram including 
the magnetic and orbital orders and the superconductiv- 
ity. For U > U' , the s±-wave superconductivity is real- 
ized near the ISDW with q ~ (tt, it) phase. On the other 
hand, for U < U', the s ++ -wave superconductivity ap- 
pears near the ferro-orbital ordered phase. The s±-wave 
pairing is mediated by the spin fluctuations, while that 
the s++-wave pairing is mediated by the orbital fluctua- 
tions. 

For U > U', the gap function for the s±-wave pairing 
changes its sign between the hole pockets and the elec- 
tron pockets and the most dominant contribution of the 
gap function is the d x 2_ y 2 orbital diagonal component. 
This is qualitatively consistent with the results based on 
the 5-band Hubbard model 25-30 . However, the d x i_ y i 
diagonal component of the gap function in our 16-band 
d-p model have much larger value than the other matrix 
elements in comparison with the results based on the 
5-band Hubbard model 26,28 . This may be because the 
outer hole Fermi surface which has mainly d yz and d zx 
orbital character obtained by the d-p model is almost cir- 
cular, but that obtained by the 5-band Hubbard model 
is diamond shape 25 . Therefore, the nesting effect which 
enhances the spin fluctuations becomes weak in our d-p 
model, and the resulting components of A(k) related to 
dyz, d zx orbitals have smaller values. 

For U < U' , the gap function in the orbital representa- 
tion for the s++-wave pairing dose not change its sign in 



k space. In other words, the on-site pairing is much larger 
than the off-site pairing in the real space. This is similar 
to the conventional phonon-mcdiated superconductivity. 
However, the gap functions have different signs between 
orbitals in contrast to the conventional phonon-mediated 
superconductivity. We have shown that this sign change 
of the gap functions between orbitals is due to the ef- 
fect of the pair transfer interaction J /53 i 54 . It is noted 
that the s++-wave state has been observed also in the 
one-dimensional 2-band Hubbard model in the same pa- 
rameter region with U < U'—. 

It seems that the both s-t-wave and the s ++ -wave 
states with full superconducting gaps are consistent with 
various experiments such as, the NMR relaxation rate, 
the Knight shift, the ARPES, the magnetic penetration 
depth measurements, although the sign of the gap func- 
tion has not been directly observed there. However, ac- 
cording to the recent theoretical studies of the nonmag- 
netic impurity effects 5 -^, Anderson's theorem is violated 
for the s±-wave superconductivity in contrast to the ex- 
perimental results of very weak T c -suppression in Fe site 
substitution^ and neutron irradiation^ 6 -. Since it can 
be considered that the impurity potential by the Fe-site 
substitution is diagonal and local in the orbital basis ac- 
cording to the first principle calculation^, it is expected 
that the s++-wave state observed in the present study 
is more robust against the nonmagnetic impurity effects 
than the s±-wave state. 

In addition to the Coulomb interaction, we have also 
discussed the effects of the coupling g between the elec- 
tron and the A\ g phonon within the RPA. It has been 
found that the s ++ -wave pairing realized in the unrealis- 
tic parameter region with U < U' for g — is enhanced 
due to the effect of g and can be expanded over the real- 
istic parameter region with U > U' for a realistic value 
of g. In the first principle calculations for iron-based su- 
perconductors in conjunction with the Migdal-Eliashberg 
theory, the electron-phonon coupling is found to be too 
small to obtain high T c observed in experiments^ 3 -. The 
effect of the Coulomb interaction, however, has not been 
discussed there. In the present study, the cooperative ef- 
fect of the Coulomb interaction and the electron-phonon 
coupling is crucial for the enhancement of the orbital 
fluctuations which induce the s++-wave superconductiv- 
ity. Recently, the large isotope effects on the transition 
temperatures for both the SDW and the superconduc- 
tivity have been observed 5 ^,. This experimental result 
implies that not only the Coulomb interaction but also 
the electron-phonon coupling plays crucial effects on the 
electronic states for iron-based superconductors. 

In early theoretical studies for the copper oxide su- 
perconductors, the effect of the Coulomb interaction be- 
tween the d and p electrons U p d was studied by several 
author s 60 ' 61 . According to the RPA study based on the 
d-p model with the single d x 2_ y 2 orbital, U pc i enhances 
the charge fluctuations with q = (0, 0) and the s-wave 
superconductivity is realized due to the effect of charge 
fluctuations 6 ^. In addition, the 1/iV-expansion ap- 
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proaches (N is the spin-orbital degeneracy) revealed that 
the strong correlation effect enhances the charge fluc- 
tuations together with the s-wave superconductivity 61 . 
Therefore, it is expected that, in the present d-p model 
with multi d orbitals, U p d enhances the charge-orbital 
fluctuations which induce the s++-wave superconductiv- 
ity. The explicit calculations based on the d-p model 
including not only the on-site Coulomb interaction but 
also the inter-site Coulomb interaction U p d together with 
the electron-phonon coupling g are now under way. 
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